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1  Introduction 


This  aspect  of  the  project  is  concerned  with  coordinating  the  sensors  in  an  underwater 
acoustic  network  in  order  to  collaboratively  track  an  acoustic  source.  Measurements  are 
taken  at  each  sensor  node,  and  in  order  to  obtain  the  best  accuracy,  the  measurements 
should  be  jointly  processed  or  fused.  This  requires  communication  and  coordination  among 
the  nodes.  At  the  same  time,  underwater  communication  is  notoriously  challenging.  Chan¬ 
nel  conditions  change  rapidly  and  high  data-rate  communications  are  generally  not  possi¬ 
ble.  Consequently,  protocols  and  mechanisms  must  be  used  which  can  adapt  to  the  time- 
varying  and  unreliable  communication  medium. 

Tracking  strategies  that  rely  on  relaying  all  measurements  to  a  central  location  become  in¬ 
feasible,  both  because  of  communication  overhead  and  robustness  concerns.  We  therefore 
focus  on  decentralized  methods  in  which  the  microprocessor  system  attached  to  each  sensor 
performs  tracking.  These  local  tracking  algorithms  share  information  between  neighbour¬ 
ing  sensors,  distilling  the  raw  measurements  into  essential  summary  statistics  so  that  the 
communication  burden  is  reduced. 

In  this  first  year  of  the  project,  we  have  focused  on  developing  distributed  algorithms  for 
single-target  tracking.  In  the  second  year,  we  will  migrate  our  attention  to  the  multi-target 
scenario.  We  have  strived  to  keep  our  research  activities  tied  to  realistic  measurement 
environments  by  repeatedly  assessing  performance  using  experimental  datasets.  The  main 
application  of  interest  involvess  bearings-only  measurements,  but  we  have  also  considered 
other  measurement  scenarios  in  order  to  make  a  more  generalizable  reasearch  contribution. 
We  have  explored  and  developed  two  algorithmic  methodologies,  concentrating  on  tracking 
solutions  that  can  adequately  address  the  non-linearities  both  in  target  dynamics  and  sensor 
measurements. 

2  Problem  formulation 


Our  goal  is  to  track  the  state  of  a  single  noise  source  using  measurements  gathered  from 
a  distributed  network  of  sensors.  We  denote  the  state  of  the  noise  source  at  time  t  by 
xf.  Typically  the  state  is  a  four-dimensional  vector,  with  two  dimensions  for  the  source 
position  and  two  dimensions  for  the  source  velocity.  The  state  evolves  according  to  a 
dynamic  model, 

Xf+l  =<?(xf)+vr ,  (1) 

where  g(-)  is  a  (possibly  non-linear)  mapping  describing  the  evolution  of  the  state  from  one 
time  step  to  the  next,  and  \t  is  a  random  vector  of  process  noise.  This  model  is  very  general 
and  can  be  used  to  capture  a  wide  range  of  noise  source  dynamics.  This  sort  of  first-order 
model  (where  xH  |  only  depends  on  xt)  is  standard  in  the  tracking  and  filtering  literature. 
Examples  of  specific  dynamic  models  are  discussed  in  further  detail  in  Section  4  below. 
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Measurements  are  obtained  by  a  network  of  sensors.  We  assume  there  are  n  sensors,  and 
the  measurement  obtained  at  sensor  j  at  time  t  is  given  by 


zr0)  =hj,t(xt)  +  Wj,t  ,  (2) 

where  hjj(-)  is  the  function  describing  the  mapping  from  the  current  noise  source  state  to 
the  observation  at  node  j  at  time  t,  and  w7.,  is  additive  measurement  noise  at  node  j  at  time 
t.  Again,  the  measurement  model  is  general.  Of  particular  interest  in  this  project  is  the 
case  where  hj_t(xt)  is  the  true  bearing  angle  from  the  sensors  position  at  time  t  to  the  noise 
source  position. 

3  Methodology 


In  our  proposed  approach,  each  node  runs  its  own  local  instance  of  the  tracking  algorithm. 
Rather  than  having  each  node  operate  in  isolation,  with  node  j  updating  its  state  estimate 
using  only  its  own  measurements,  our  aim  is  develop  schemes  under  which  sensors  share 
information  about  their  local  measurements  so  that  the  tracking  performance  of  the  overall 
system  is  superior  to  that  of  any  individual  sensor.  The  main  issues  to  be  addressed  are 
then: 


•  What  information  should  be  communicated? 

•  What  tracking  algorithm  should  be  used  at  each  node? 

Below  we  describe  our  approach  to  these  issues.  With  each  node  running  a  local  tracking 
algorithm,  communication  plays  the  role  of  synchronizing  the  state  estimates  across  all 
nodes  in  the  network  so  that,  ideally,  all  nodes  have  the  same  estimate  that  would  be  com¬ 
puted  by  a  single  tracking  algorithm  that  had  direct  access  to  all  of  the  measurements.  We 
use  gossip  algorithms  (discussed  in  Section  3.1)  to  diffuse  information  across  the  network 
and  drive  the  state  estimates  at  each  node  to  a  consensus. 

For  the  local  tracking  algorithm  used  at  each  node,  we  investigate  two  options:  the  particle 
filter  and  the  ensemble  Kalman  filter  (discussed  in  Sections  3.2  and  3.3).  Particle  filters  are 
attractive  when  the  dynamic  model  g(-)  is  non-linear  (e.g.,  if  the  noise  source  may  be  highly 
manoeuvring),  and  when  the  measurement  model  hjj(-)  is  non-linear,  as  is  certainly  the 
case  for  the  bearings  measurement  model  considered  in  this  project.  However,  the  overhead 
associated  with  communicating  the  particle  representation  is  non-trivial.  Variations  of  the 
ensemble  Kalman  filter  make  it  possible  to  reduce  the  communication  overhead,  possibly 
at  the  cost  of  some  performance  degradation. 
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3.1  Gossip  algorithms  for  distributed  averaging 

We  begin  by  discussing  gossip  algorithms,  which  we  use  to  synchronize  and  spread  infor¬ 
mation  across  the  network.  Gossip  algorithms  are  a  class  of  distributed  message -passing 
algorithms  for  aggregating  data  over  a  network  [1,  2,  3];  see  [4]  for  a  survey.  We  focus  on 
gossip  algorithms  which  solve  the  distributed  averaging  problem :  each  node  j  =  1 .....  n 
initially  has  a  value  yj( 0),  and  the  goal  of  the  network  is  for  all  nodes  to  obtain  the  average 
of  the  initial  values,  y  =  ^  Y!)=\  >y(0).  In  order  to  reach  this  goal,  the  nodes  repeatedly 
transmit  messages  and  modify  their  local  estimate  of  the  average  based  on  messages  they 
receive.  This  process  can  occur  synchronously  (all  nodes  transmit,  and  then  they  update 
simultaneously  before  the  next  round  of  communication)  or  asynchronously  (where  nodes 
update  whenever  they  receive  a  message  and  then  broadcast  their  new  result  periodically). 

In  the  most  general  case,  after  k  updates  have  occurred  across  the  network,  each  node  has  a 
current  estimate  yj{k)  for  the  average  value.  At  the  next  update/iteration,  suppose  that  node 
i  receives  messages  from  a  subset  of  nodes  Nj(k)  C  (1,2, ...  ,n}.  For  each  node  j  e  Nj(k), 
the  message  received  by  node  i  from  j  contains  the  value  }'j(k).  Then  node  i  updates  its 
new  value  to 

yi(k+l)  =Aij{k)yi{k)+  £  Aij(k)yj(k) ,  (3) 

jeNi(k) 

where  the  weights  (Ay(k)}”=1  typically  satisfy  (i)  ALj(k)  >  0;  (ii)  Aa(k)  >  0;  and  (iii) 
Y!j=\  Aij(k)  —  1.  These  conditions  effectively  state  that  (a)  nodes  form  their  new  estimate 
by  taking  a  weighted  average  of  their  old  estimate  with  the  estimates  received  from  their 
neighbours,  (b)  they  always  give  some  weight  to  their  own  estimate,  and  (c)  the  total  weight 
given  to  the  estimates  forming  v/(k+  1)  is  equal  to  1. 

Of  particular  interest  in  this  project  is  the  case  where  each  node  is  only  ever  able  to  com¬ 
municate  directly  with  a  subset  of  nodes,  and  this  connectivity  is  represented  as  a  graph 
G  =  (V,E)  with  vertices  V  =  {l,...,n}  corresponding  to  the  nodes,  and  with  an  edge 
( i ,  j)  e  E  if  node  i  can  receive  a  message  from  node  j.  In  order  for  there  to  be  any  hope  of 
the  entire  network  converging  to  a  consensus  on  the  average,  the  graph  of  possible  connec¬ 
tions,  G,  should  be  strongly  connected :  there  should  be  a  path  (possibly  involving  multiple 
hops)  from  every  node  i  to  all  other  nodes  j  e  V  \  {/}  in  the  network. 

For  most  gossip  algorithms  it  can  be  guaranteed  that  all  nodes  converge,  asymptotically,  to 
a  consensus  on  the  average,  and  the  rate  of  convergence  depends  on  structural  properties 
of  G.  More  specifically,  after  k  iterations  we  have  (with  high  probability,  in  the  case  of 
asynchronous  randomized  gossip), 

t  (y,(k)  -  yf  <  Pk  t  MO) -y)2,  (4) 

i=l  i=l 

where  0  <  p  <  1  is  a  constant  related  to  the  connectivity  of  the  graph  G.  Roughly  speaking, 
if  G  is  very  well  connected  (e.g.,  there  are  many  short  paths  between  every  pair  of  nodes), 
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then  p  is  closer  to  0  and  the  algorithm  converges  faster.  Similarly,  if  there  is  a  probability 
that  packets  are  not  successfully  received  then  this  increases  p  (bringing  it  closer  to  1, 
corresponding  to  slower  convergence). 

Below  we  will  discuss  how  gossip  algorithms  can  be  used  as  a  subroutine  within  distributed 
tracking  algorithms  to  synchronize  the  state  estimates  at  each  node.  For  distributed  track¬ 
ing,  each  node  has  its  own  local  representation  of  the  state  estimate,  and  we  would  like 
to  keep  these  representations  synchronized.  Nodes  will  update  these  local  representations 
using  their  local  measurement,  and  then  they  will  perform  some  gossip  iterations  in  order 
to  re- synchronize  (equivalently,  to  spread  the  information  each  node  learned  from  its  local 
measurement  across  the  network).  If  the  state  estimates  are  well- synchronized  to  begin, 
and  if  the  update  using  local  measurement  does  not  change  the  state  estimate  too  much, 
then  after  a  few  gossip  iterations  they  will  be  synchronized  again. 


3.2  Distributed  particle  filtering 


The  dynamic  model  (1)  defines  the  distribution  p(xt  |xf_i)  of  the  state  at  time  t  given  the 
state  at  time  t  —  1.  Similarly,  the  measurement  model  (2)  defines  the  likelihood  p(zt\xt)  of 
an  observation  zt  given  that  the  state  is  xf.  Given  these  models,  the  Bayes-optimal  recursion 
for  the  posterior  state  estimate  is  given  by  the  equations: 


Predict: 

Update: 


p{xt  |z1:f_ i)  = 
p(xt\zht)  = 


Xf  |Xf_!  )  p(X;_i  |z1:f_!  )  d,Xf—\ 

p(zt\xt)  p(xt\zUt_i) 
p(zr|zi;/_i) 


(5) 

(6) 


where  Zj _  i  denotes  the  collection  of  all  measurements  from  times  1  up  to  /  —  1 .  Except  for 
under  very  special  circumstances  (e.g,.  unless  the  models  g(-)  and  h(-)  are  both  linear  and 
the  noise  terms  \t  and  w,  are  Gaussian),  these  recursions  cannot  be  computed  analytically. 
Hence,  in  most  cases,  we  must  resort  to  some  approximation. 


Particle  filters  approximate  the  posterior  distribution  p(xt \z\:t)  using  a  set  of  m  weighted 
particles,  {xf\  \  specifically, 

An(yr|zi:r)  =  £  W,(0  8{xt  -  xj0)  ,  (7) 

i=  1 

and  where  <5(-)  is  the  Dirac  delta  function.  Each  particle  xW  can  be  viewed  as  a  candidate 
hypothesis  for  the  current  state  value,  and  the  corresponding  importance  weight  w\1'1  is  an 
approximation  to  the  posterior  density  at  x\l\  The  particles  are  updated  from  one  time  step 
to  the  next  (i.e.,  from  xp  to  xj" , )  by  simulating  the  noise  source  dynamic  model  (1),  and 

then  the  new  importance  weights  w\'lx  are  updated  based  on  the  new  position  x|^  and  the 
new  measurements  zt+\.  For  more  background  on  particle  filters,  see  [5,  6]. 
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For  distributed  particle  filtering,  we  consider  the  case  where  each  node  maintains  its  own  set 
of  weighted  particles.  Let  zf  =  {ztJ'  :  j  =  1, ...  ,n}  denote  the  collection  of  measurements 
from  all  nodes  at  time  t.  We  assume  that  measurements  at  different  nodes  are  conditionally 
independent  given  the  state.  The  joint  likelihood  then  factorizes  into 

P{ z,v|x)  =  f\p{z(tj)\xt)  • 

7=1 

Then,  to  incorporate  information  from  all  of  the  measurements,  the  weight  of  the  ith  particle 
should  be  set  to 

w?  oc  exp  (  logp(z}j) |xfW)  j  ,  (8) 

where  the  proportionality  constant  only  depends  on  the  particle  x,l\  and  not  on  the  mea¬ 
surements.  Hence,  in  the  distributed  setting,  to  fuse  measurements  gathered  at  all  sensors 
we  need  to  compute  the  sum  in  (8).  Of  course,  after  rescaling  each  term  in  the  sum  by  n, 
this  is  equivalent  to  computing  an  average, 

i^n]ogp(z[j}\xll))  , 

and  so  we  can  use  the  gossip  algorithms  described  in  Section  3.1  with  the  value  at  node  j 
initialized  to  =  n  log p(z!t  n  x)^). 

The  performance  of  particle  filters,  in  general,  is  related  to  the  number  of  particles  used.  As 
the  number  of  particles  m  tends  to  infinity,  the  particle  approximation  (7)  converges  to  the 
posterior.  On  the  other  hand,  in  the  setting  of  distributed  particle  filtering,  we  need  to  reach 
a  consensus  on  each  particle  weight.  Thus,  each  gossip  message  is  of  dimension  m.  Conse¬ 
quently,  we  have  a  tradeoff  in  distributed  particle  filtering,  where  increasing  the  number  of 
particles  m  should  improve  the  tracking  accuracy  at  the  cost  of  increased  communications. 
Increasing  m  also  increases  the  amount  of  computation  each  node  must  perform  in  order  to 
update  the  particle  locations  and  weights,  but  we  believe  this  is  a  secondary  concern  in  the 
context  of  this  project,  since  the  underwater  communication  medium  is  the  main  bottleneck 
limiting  the  number  of  particles  that  may  be  practical. 

3.3  Distributed  ensemble  Kalman  filtering 

The  distributed  ensemble  Kalman  filter  is  an  alternative  approach  to  tracking,  which  we 
consider  promising  because  it  allows  us  to  reduce  the  communication  overhead  [7].  The 
(centralized)  ensemble  Kalman  filter  [8]  uses  a  set  of  samples,  which  we  denote  by  xj‘\ 
to  approximate  the  first  two  moments  of  the  state  density.  Two  particular  versions  of  the 
ensemble  Kalman  filter — the  ensemble  square  root  filter  [9]  and  the  deterministic  ensemble 
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Kalman  filter  [9] — are  particularly  attractive  for  use  in  the  distributed  setting  because  their 
updates  dramatically  reduce  the  amount  of  communication.  In  these  filters  each  of  the 
samples  is  updated  using  a  Kalman  filter-like  recursion, 

Forecast:  x^  =  g(x^j)  (9) 

Update:  x^“  =  xf  +  K(y  —  Htxf)  +  T(x^  —  x{)  ,  (10) 

where  x{  =  ^E^Li  k(x/-i)  's  the  average  of  the  forecasted  samples  (i.e.,  the  average  state 
after  passing  each  sample  through  the  dynamic  model),  H,  is  the  linear  measurement  ma¬ 
trix.  K  =  p/h/  (H/P/h/  +  R)  1  represents  the  Kalman  gain,  in  which  R  is  the  measure¬ 
ment  noise  covariance  matrix  and  P/  =  ^-j-  Ef- 1  (xjM,/  —  x/)(x|,!^  —  x/)7  is  the  prediction 
covariance  matrix,  estimated  from  the  samples.  The  operator  T  =  (I  —  KHr)1//2  for  the 
ensemble  square  root  filter,  where  I  is  the  identity  matrix,  and  T  =  I  —  ^KH;  for  the  de¬ 
terministic  ensemble  Kalman  filter.  These  equations  are  suitable  for  linear  measurement 
models,  but  for  non-linear  measurement  models,  we  perform  a  linearization  of  each  sen¬ 
sor’s  measurement  function  hj  t(x W)  to  derive  a  matrix  Hf  (due  to  the  linearization,  this 
matrix  can  vary  for  each  sample  ;). 

The  main  challenge  in  formulating  a  distributed  update  step  lies  in  the  fact  that  the  forecast 
and  update  equations  (9)-(10)  require  the  knowledge  of  the  complete  set  of  measurements, 
measurement  functions  and  measurement  noise  statistics.  Sensor  nodes  only  have  access  to 
their  local  measurements  and  may  only  be  aware  of  their  own  measurement  model  parame¬ 
ters.  In  [7],  we  addressed  this  challenge  by  expressing  the  update  equation  in  an  alternative 
information  filter  form.  We  identify  summary  statistics  that  can  be  computed  by  global 
network  averaging  of  locally-computable  information  using  gossip  procedures. 


4  Results 

4.1  Distributed  particle  filtering 

Next  we  summarize  the  results  obtained  for  distributed  particle  filtering.  A  preliminary 
version  of  our  results  are  presented  in  the  conference  paper  [10],  and  an  extended  version 
is  in  preparation  to  be  submitted  to  a  journal  shortly. 

To  evaluate  the  performance  of  the  proposed  gossip-based  distributed  particle  filtering  so¬ 
lution,  we  use  data  from  an  at-sea  experiment  conducted  in  October  2012  in  the  Emerald 
Basin  area  of  the  Scotian  Shelf  [1 1].  The  experimental  setup  involves  five  acoustic  sensors, 
a  single  acoustic  source,  and  an  echo-repeater  in  the  water.  We  focus  on  localizing  and 
tracking  the  source  using  measurements  obtained  at  the  sensors.  Figure  1  depicts  the  tra¬ 
jectory  of  the  source  along  with  the  trajectories  of  the  sensors,  which  were  freely  drifting. 
The  experiment  lasts  for  roughly  180  minutes,  with  each  sensor  providing  measurements 
once  per  minute.  During  this  time,  the  noise  source  makes  two  cycles. 
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Figure  1:  (a)  Trajectories  of  the  source  and  the  sensors  in  the  experiment  considered,  (b) 
Time-averaged  squared  error  vs.  reception  probability  p  for  a  particle  filter  using  both  bear¬ 
ings  and  range  measurements,  (c)  Time-averaged  disagreement  vs.  reception  probability  p. 
The  results  shown  correspond  to  1 000  Monte  Carlo  trials. 


We  consider  two  dynamic  models.  The  first  models  the  noise  source  motion  at  each  time- 
step  as  either  remaining  true  with  a  constant  velocity  (with  probability  pcv)  or  making 
a  turn  with  constant  turning  rate  (with  probability  1  —  pcv)-  The  approach  combines  a 
constant  velocity  model  (where  the  source  moves  along  a  straight  line)  and  a  coordinated 
turn,  so  we  refer  to  it  as  the  CV-CT  model.  The  second  model  we  consider  is  a  simple 
random  walk  model  (which  we  refer  to  as  the  RW  model),  where  the  motion  of  the  noise 
source  at  each  time  step  is  random  and  isotropic  about  its  previous  state.1 

Our  preliminary  investigation,  reported  in  [10],  considered  the  case  where  each  sensor  ob¬ 
tains  noisy  bearings  measurements  as  well  as  a  noisy  range  measurement.  For  this  setting, 

'in  both  dynamic  models,  we  model  the  process  noise  and  measurement  noise  w,  as  being  zero-mean 
Gaussian,  with  a  diagonal  covariance  matrix  whose  parameters  were  estimated  from  the  data  [10].  In  the 
CV-CT  model,  we  obtain  estimates  for  the  turn  probability  and  the  turn  radius  from  the  dataset. 
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we  simulate  the  case  where  each  time  a  sensor  transmits  a  message  it  is  received  at  each 
other  sensor  independently  with  some  probability  p,  and  we  study  the  tracking  performance 
as  a  function  of  the  reception  probability.  This  is  shown  in  Figure  1  (b),  for  a  filter  with  1600 
particles,  and  Figure  1(c)  shows  the  time-averaged  disagreement2  between  nodes  as  a  func¬ 
tion  of  p,  under  the  CV-CT  model.  Figure  1  shows  that  the  message  reception  probability 
has  a  non-trivial  impact  on  the  performance.  The  time-average  squared  error  clearly  in¬ 
creases  when  the  reception  probability  is  lower.  As  the  probability  of  receiving  messages 
from  other  sensors  approaches  zero,  the  performance  corresponds  to  that  of  each  sensor 
performing  tracking  in  isolation.  Clearly  collaboration  is  beneficial,  both  for  reducing  the 
mean  time-averaged  square  error  as  well  as  reducing  the  error  variance.  From  Figure  1(c), 
we  also  see  that  as  there  is  less  collaboration  (i.e.,  when  sensors  receive  fewer  messages 
from  their  neighbours)  there  is  greater  disagreement  among  the  state  estimates  at  different 
sensors. 

Next  we  focus  on  the  performance  when  sensors  only  use  bearings  measurements.  Fig¬ 
ure  2(a)  shows  the  root  mean-squared  error  (RMSE)  as  a  function  of  the  number  of  parti¬ 
cles  used  in  a  centralized  particle  filter  (i.e.,  one  having  access  to  the  measurements  from  all 
sensors),  and  Figure  2(b)  shows  the  computation  time  per  filter  update  as  a  function  of  the 
number  of  particles.  Note  that  the  horizontal  axis  of  these  figures  is  linear  over  the  range 
from  10  to  100,  and  then  it  increases  geometrically  from  100  to  1600.  The  experiment 
was  conducted  on  a  computer  with  a  1.7GHz  Intel  i7  processor  with  8GB  of  RAM.  The 
performance  of  the  centralized  filter  gives  a  benchmark  for  what  we  can  hope  to  achieve 
in  the  distributed  setting.  From  Figure  2(b)  we  observe  that  the  computation  time  scales 
linearly  with  the  number  of  particles  used.  Figure  2(c)  shows  the  RMSE  over  time  of  the 
centralized  particle  filter  as  a  function  of  time  for  the  two  different  dynamic  models.  The 
average  error  is  lower  for  the  CV+CT  model,  which  makes  stronger  assumptions  about  the 
noise  source  dynamics,  but  the  CV+CT  model  does  have  higher  peak  error. 

To  better  understand  the  cause  of  these  peak  errors  under  the  CV+CT  model,  we  observe 
the  timing  of  these  errors  coincides  with  times  when  the  noise  source  is  making  a  turn 
towards  the  sensors,  as  shown  in  Figure  3.  During  these  time  periods,  the  noise  source 
trajectory  is  roughly  co-linear  with  the  sensor  positions.  Bearings-only  tracking  is  known 
to  be  especially  challenging  when  the  source  moves  in  a  line  which  is  co-linear  with  the 
sensors,  since  then  the  bearings-measurements  do  not  change  significantly  although  the 
source  may  have  moved  a  large  distance.  This  suggests  that  it  is  also  important  to  consider 
the  problem  of  sensor  placement. 

2  Here  the  disagreement  is  the  average,  over  all  pair  of  nodes,  of  the  square  difference  between  the  nodes’ 
maximum  a  posteriori  state  estimate. 
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Figure  2:  Performance  of  the  centralized  particle  filter  with  the  number  of  particles  rang¬ 
ing  from  10  to  1600.  (a)  Root  mean-squared  error  vs.  number  of  particles,  (b)  Runtime 
vs.  number  of  particles,  (c)  Example  tracking  performance  under  the  two  different  motion 
models  for  the  centralized  particle  filter  with  1600  particles. 


4.2  Distributed  ensemble  Kalman  filtering 

In  order  to  examine  the  performance  of  the  proposed  distributed  ensemble  Kalman  filter, 
we  simulate  a  nonlinear  radio -frequency  (RF)  tomography  tracking  system  (see  [12]  for  a 
complete  description  of  the  RF  tomography  model).  Twenty  four  nodes  are  deployed  at  the 
boundary  of  a  square  region  of  50m  x  50m.  The  sensor  nodes  track  a  single  target  for  50 
seconds,  making  a  set  of  measurements  every  second.  Sensor  nodes  can  communicate  with 
each  other  if  their  separation  is  less  than  a  transmission  range  of  15m.  The  target  motion  is 
modeled  by  the  nearly  coordinated  turn  model,  which  assumes  unknown  cartesian  velocity 
but  known  turn  rate.3  Table  1  compares  the  tracking  accuracy  of  the  distributed  ensemble 
Kalman  filters  we  provided  and  state-of-the-art  non-linear  distributed  tracking  alternatives. 

3We  use  the  RF  tomography  measurement  model  from  [12],  with  parameters  oy  =  0.05  and  (j)  =  5,  and 
measurement  noise  standard  deviation  of  0.50.  The  standard  deviation  of  the  target  motion  is  0.25. 
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Figure  3:  RMSE  of  the  particle  filter  under  the  CV+CT  model.  Red  segments  indicate 
times  where  the  RMSE  peaks.  These  coincide  with  times  when  the  noise  source  turns 
towards  the  sensors  and  its  trajectory  is  roughly  co-linear  with  the  sensors. 


As  we  reduce  the  number  of  scalars  that  each  node  is  permitted  to  transmit,  we  see  that  the 
proposed  algorithms  clearly  outperform  the  alternative  methods. 

5  Summary  and  ongoing  work 


We  have  investigated  the  problem  of  distributed  target  tracking.  Our  main  application 
of  interest  involves  bearings-only  measurements  using  a  network  of  multiple  underwater 
acoustic  sensors.  The  challenges  of  underwater  communication  have  motivated  us  to  focus 
on  distributed  tracking  algorithms  that  are  capable  of  addressing  the  non-linearities  in  both 
target  dynamics  and  sensor  measurements.  We  have  therefore  concentrated  our  attention 
on  distributed  particle  filtering  and  distributed  ensemble  Kalman  filtering. 

We  have  conducted  a  thorough  performance  assessment  of  candidate  distributed  particle 
filtering  algorithms  using  experimental  datasets.  To  the  best  of  our  knowledge,  this  is  the 
first  time  that  decentralized  particle  filters  have  been  carefully  compared  using  experimen¬ 
tal  data.  Our  research  has  highlighted  the  need  for  improving  the  distribution  of  the  sensors 
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Table  1:  Tracking  Performance  for  the  Distributed  Ensemble  Kalman  Filter 


Ave.  RMSE  ±  Std.  Dev.  (%  track  loss) 

Scalars 

1000 

500 

200 

Gaussian  Approximation  [13] 

0.17  ±0.01 

0.20  ±  0.02(2) 

0.35  ±0.18(22) 

Likelihood  Consensus  [14] 

0.33  ±  0.22 

- 

- 

Set  Membership  [15] 

0.26  ±0.18(3) 

0.33  ±0.19(7) 

0.36  ±0.25(16) 

Top-m  Selective  Gossip  [16] 

0.27  ±  0.22(2) 

0.34  ±0.20(11) 

0.70  ±  0.41(50) 

Distributed  UKF  [17] 

0.21  ±0.09(18) 

0.61  ±  0.49(74) 

0.40  ±  0.26(95) 

Distributed  ESRF 

0.19  ±0.05 

0.18  ±  0.05 

0.22  ±  0.10(6) 

Distributed  DEnKF 

0.23  ±  0.09 

0.21  ±  0.08 

0.22  ±  0.09(1) 

Centr.  PF 

0.10  ±0.01 

0.10  ±0.01 

0.10  ±0.01 

to  enhance  observability  of  the  targets,  and  the  need  to  develop  algorithms  that  can  maintain 
a  reasonable  accuracy  but  have  much  less  communication  overhead.  With  this  in  mind,  we 
have  proposed  a  novel  distributed  filtering  algorithm,  developing  techniques  for  using  gos¬ 
sip  algorithms  to  implement  ensemble  Kalman  filtering  calculations.  Compared  to  other 
state-of-the-art  non-linear  tracking  algorithms,  simulation  results  suggest  that  the  proposed 
filters  can  maintain  reasonable  accuracy  while  reducing  the  communication  requirements. 

In  the  next  year  of  the  project,  we  will  migrate  our  attention  to  the  multi-target  tracking 
problem,  which  poses  additional  data  association  challenges  and  involves  tracking  in  a 
much  higher  dimensional  state-space.  We  are  currently  developing  techniques  for  imple¬ 
menting  distributed  multi-sensor  random  finite-set  filters. 
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